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Abstract 

Using the large Nq approximation we have constructed the most general chiral reso- 
nance Lagrangian in the odd-intrinsic parity sector that can generate low energy chiral 
constants up to 0(p e ). Integrating out the resonance fields these 0(p 6 ) constants are 
expressed in terms of resonance couplings and masses. The role of rf is discussed and its 
contribution is explicitly factorized. Using the resonance basis we have also calculated two 
QCD Green functions of currents: (VVP) and (V AS) and found, imposing high energy 
constraints, additional relations for resonance couplings. We have studied several phe- 
nomenological implications based on these correlators from which let us mention here our 
prediction for the 7r°-pole contribution to the muon g — 2 factor: = 65.8(1.2) x 1CU 11 . 
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1 Introduction 

As is well known, there are two regimes where the QCD dynamics of the current correlators 
is well understood. The first one corresponds to the high energies where the asymptotic 
freedom allows to use the perturbative approach in terms of the strong coupling constant 
a s and where the asymptotics of the correlators for large euclidean momenta is governed by 
operator product expansion (OPE). The second well understood region is that of low external 
momenta where the dynamics is constrained by the spontaneously broken chiral symmetry. 
As a consequence, the dominant contributions to the correlators and related amplitudes of 
the processes under interest come from the octet of the lightest pseudoscalar mesons (tt, K, 
j]) which are the corresponding (pseudo)Goldstone bosons (GB). The correlators can be 
studied here by means of Chiral Perturbation Theory (ChPT)[Q [3], which is the effective 
Lagrangian field theory for this region, in terms of systematic simultaneous expansion in 
powers (and logs) of the momenta and quark masses. The applicability of ChPT extends up 
to the hadronic scale Ah ~ lGeV which corresponds to the onset of non-Goldstone resonances 
and where the ChPT expansion fails to converge. 
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OPE and ChPT provides us with asymptotic behaviour of the correlators in different 
regimes, however, both these approaches need further non-perturbative long-distance piece 
of information which is not known from the first principles, namely the values of the vacuum 
condensates for OPE and the values of the effective low-energy constants (LECs) for ChPT. In 
the latter case the LECs parameterize our lack of detailed information on the non-perturbative 
dynamics of the degrees of freedom above the hadronic scale Ah and are connected with the 
order parameters of the spontaneously broken chiral symmetry. The predictivity of ChPT 
heavily relies on their determination. At the order 0(p e ), which corresponds to the recent 
accuracy of the NNLO ChPT calculation (for a comprehensive review and further references 
see [1]), 90+4 LECs in the even intrinsic parity sector [SJ |6] and 23 LECs in the odd sectoi0 
|8] appear in the effective Lagrangian. Though only special linear combinations of them 
are relevant for particular physical amplitudes, the uncertainty in their estimation is usually 
the weakest point of the interconnection between the theory and experiment. 

Dispersion representation of those correlators which are order parameters of the chiral 
symmetry breaking (and therefore do not get any genuine perturbative contribution) enables 
to make use of information on the asymptotics both in the low and high energy regions and to 
relate the unknown LECs to the properties of the corresponding spectral functions in terms 
of the chiral sum rules [21 El [101 [III 02] . These are usually assumed to be saturated by the 
low-lying resonant states; such an assumption (known as resonance saturation hypothesis) 
connects the LECs to the phenomenology of resonances in the intermediate energy region 
lGeV < E < 2GeV. Though the inclusion of only finite number of resonances has been 
questioned in the literature [El [14], it proved to be consistent in the 0(p 4 ) case with other 
phenomenological determinations of LECs. 

The necessary ingredient of the resonance saturation approach to the determination / esti- 
mation of LECs is the phenomenological information on the physics of the lowest resonances. 
It can be conveniently parameterized by means of suitable phenomenological Lagrangian. 
Along with the chiral symmetry the guiding theoretical principles for its construction are 
those based on the large Nc expansion of QCD [15J. Within the leading order in 1/Nc the 
correlators of the quark bilinears are given by an infinite sum of contributions of narrow 
meson resonance states the mass of which scales as 0(N C ) and the interaction of which 
is suppressed by an appropriate power of 1/^/Nc- Such a large Nc representation of the 
correlators can be reconstructed using effective Lagrangian including GB and infinite 
tower of resonance fields with couplings of the order 0(N C n ^ 2 ) according to the number n 
of the resonance fields in the interaction vertices. The 1/Nc expansion is equivalent to the 
quasi-classical expansion, thus at the leading order only the tree graphs contribute and each 
additional loop is suppressed by one power of 1/Nc- 

Though Coo is not known from the first principles, the information on the large Nc 
hierarchy of the individual operators together with general symmetry assumptions allows one 
to construct all the relevant terms necessary to determine the LECs in the leading order 
of the large Nc expansion up to given chiral order. The large Nc approximation of LECs 
can be then formally achieved by means of the integrating out the resonance fields from the 
Lagrangian C^. Formally one gets LECs expressed in terms of the (from the first principles 
unknown) masses and couplings of the infinite tower of resonances. 

The large Nc inspired phenomenological Lagrangian suitable for the resonance saturation 

1 These numbers of LECs are relevant for SU(3) variant of ChPT. In the 51/(2) case we get 53+4 LECs in 
the intrinsic parity even sector and 5(13) LECs in the odd sector. 
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program for LECs can be then obtained as an approximation to where only finite number 
of resonances is kept. Such a truncation of seems to be legitimate at low energies where 
the contribution of the higher resonances is expected to be suppressed. However, the lack 
of effective cut-off scale which could play here a role analogous to Ah for ChPT prevents 
us to interpret the resonance phenomenological Lagrangian as a well defined effective theory 
in the usual sense. It is rather a QCD inspired phenomenological model which should share 
as much common features with QCD as possible. The latter principle generally puts various 
constraints on its effective couplings. For instance, the finite number of resonances involved 
generally corrupts the asymptotic behaviour of the correlators required by perturbative QCD 
and OPE. However, it is natural to expect that for the correlators which are order parameters 
of the spontaneous chiral symmetry breaking that the latter behaviour extends down to the 
region of applicability of the phenomenological Lagrangian and thus it is desirable to ensure 
the correct asymptotics by means of adjusting its couplings. This is however not enough to 
fix all of them (often it is even not possible to satisfy all the OPE requirements at once by a 
finite set of resonances) therefore further phenomenological input is needed. 

At the leading order in 1/Nq, the above strategy for determination of LECs is essentially 
equivalent to the similar approach known as Minimal Hadronic Ansatz (MHA)[l6J. Within 
this approach the correlators are approximated by meromorphic function with correct pole 
structure corresponding to the resonance poles and the free parameters are fixed both by 
OPE constraints and experimental inputs. Only minimal number of resonances is taken 
into account, just those necessary to satisfy all the relevant OPE (when only the lowest 
resonances in each channel are included, the method is called Lowest Meson Dominance 
(LMD) ansatz [16] I12j . but in this case not all OPE constraints are guaranteed to be met 
\12\ I17j). Matching this ansatz to the low energy ChPT expansion enables to determine 
relevant linear combinations of LECs. 

The method based on the resonance Lagrangian is however little bit more general than 
MHA or LMD. On one hand it enables to determine (at least in principle) the individual 
LECs, not only their linear combinations connected with particular correlators, on the other 
hand it provides a natural framework for going beyond the leading order in 1 /Nq by means of 
integrating out the resonances at one loop level [18] PT9] [20l [21} [22] which also takes correctly 
into account the renormalization scale dependence of the LECs. 

The above principles of construction of phenomenological Lagrangian with resonances are 
known since 1989 when the seminal paper |23| on what is now known as Resonance Chiral 
Theory (R%T) was published. In this paper the resonance saturation of the 0(p 4 ) LECs 
has been studied systematically while the 0(p e ) LECs of the even intrinsic parity sector of 
ChPT has been systematically analyzed 17 years later in [21]. For a recent review and further 
references see [25] . 

The study of the odd intrinsic parity sector of R%T with vector resonances and corre- 
sponding saturation of the LECs for the 0(p e ) anomaly sector of ChPT started in [26J and 
|27[ [T2| 128] . where also axial vector resonances has been included and where the particular 
operator basis of the R%T Lagrangian contributing to the correlators under interest has been 
constructed. The influence of pseudoscalar resonances on the odd intrinsic parity LECs has 
been studied in [11] [T2] and corresponding part of R%T Lagrangian has been constructed in 
[29] (see also |30j). In this paper we resume this effort and construct the most general odd 
intrinsic parity sector of the R^T Lagrangian including the lowest multiplets of the vector 
V(l ), axial-vector A(l ++ ), scalar 5(0 ++ ) and pseudoscalar P(0 h ) resonances. In the 
*~ channel we introduce thus beside the GB also the lowest non-GB resonance multiplet 
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and therefore we go beyond the LMD approximation (our correlators then correspond to what 
is called in |12j as LMD+P ansatz). The resulting Lagrangian is then used for the lowest reso- 
nance saturation of the 0(p e ) anomaly sector of ChPT. We also illustrate the general strategy 
of matching the correlators with OPE on the concrete example of (VVP) and (VAS) three 
point functions and discuss related phenomenological applications. 

The paper is organized as follows. In Sect. 2 we fix our notation and remind briefly 
the principles of the construction of the Lagrangian of the RxT- Sec. 3 is devoted to the 
presentation of the complete basis of the odd intrinsic parity sector of RxT- In Sect. 4 
we discuss related phenomenological applications and in Sect. 5 we give the result of the 
resonance saturation of the odd intrinsic parity 0(p G ) LECs. A brief summary is given in 
Sec 6. The large Nq counting of the relevant operators is discussed in Appendix A and the 
operator redefinitions and reduction of the Lagrangian is studied in Appendix B. 



2 The Resonance Chiral Theory 

In what follows we will work in the chiral limit. The standard basic building block which 
includes the octet of GB (here we assume that rj has been already integrated out from our 
effective Lagrangian, for details see Appendix [A]) is 



u{4>) = exp [i-^,) » O 1 ) 
where (j) = "^^ a< / ,c \ ^ being a standard Gell-Mann matrix and 

m = * -> 0+ 76 % K ° • ( 2 ) 
V *° ->/ 

One can form the basic covariant tensors [31], [5] 

u^ = u^ = i{u)(d^ - ir^)u - u (<9 M - il^vt} , 
X± = u 1 X u 1 ± u x t u , 

h^u = Vfj,u v + VuUfj, , (3) 

with x = 2So(s + ip), where s and p stand for the scalar and pseudo-scalar external sources. 
Vector source and axial-vector source are then related to the right and left sources 
r M and by relations = ^(r M + i^) and a M = ^(r M — B 1 ) respectively, and F£ U R the 
corresponding left and right field-strength tensors: 

Fg = d»r v - d u r>* - i[r», r u ] , Ff = dH v - d v l» - l u ] . (4) 

The covariant derivative is defined by 

V^X = d fl + [T fl ,X], (5) 

where the chiral connection is 

r ^ = \ W{dv ~ ir v) u + u (. d ft ~ ■ ( 6 ) 
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Inspired by the large Nc limit the GB couple to massive U(3) multiplets of the type 
V(l ), A(l ++ ), S(0 ++ ) and P(0 h ), denoted generically as a nonet field R. This field can 
be decomposed into octet Rg and singlet Ro via 

The explicit form of the vector multiplet V(l ) is 



^1/ 



(and similarly for other types). We use here the antisymmetric tensor field for description of 
the spin-1 resonances. The reason is that though it is in principle equivalent to the Proca field 
formalism (see [32] and [33J for the general discussion of the equivalence at the order 0(p 4 ) 
and 0(p®) respectively and [34J for particular discussion of the one- loop equivalence), the 
antisymmetric tensor field naturally couples to the lowest order 0(p 2 ) chiral building blocks 
without derivatives and therefore it does not require additional contact terms necessary to 
compensate the wrong high energy behaviour of amplitudes and form factors under interest. 
Moreover, when using the Proca field without such contact terms it is not possible to saturate 
the 0(p 4 ) LECs in the even intrinsic parity sector and for the similar reason also the LECs 
in the 0(p e ) odd intrinsic parity sector. 

According to the large Nc counting of interaction vertices with resonances we can organize 
the Lagrangian £r x t of RxT as an expansion in the number of resonance fields, 

&RXT = &GB + £-RR,km + C,R + CrR> + C-RR'R" + ■■■ (9) 

Here Cqb contains only Goldstone bosons and external sources and includes terms with the 
same structure as the usual SU(3)l x SU(3)r ChPT Lagrangian, but the coupling constants 
are generally different. The resonance kinetic terms Cnn^n, which are of the order 0(N C ), 
have the form 

CrrM* = -\{^R^aR av ) + jM 2 R (R^R^) + l(V a R , V a R'} -\M 2 RI {R!R!) , (10) 

where R stands for V^ p and A ftv while R' stands for S and P. The terms Cr, Crri and 

£rr'R" collect the interaction vertices linear, quadratic and cubic in the resonance fields, 
respectively. 

There is also another type of expansion for Cr x t- It is based on the ordering according 
to the contribution to chiral coupling constants. Within this counting, the resonance fields 
are effectively of the order 

R = 0( P 2 ), (11) 

while the chiral building blocks with GB only are counted in a usual way. For Cqb it is 
therefore just the usual chiral power counting. Combining this with the large Nc expansion 
dH) we can write 

r n .-/-( 2 ' -ir< 4 ' i r( 6 ) , r W , r (6) r (6) r (6) ~(6) h9 , 

>~RxT — >-GB GB L 'RR,kin L 'RR,kin R '-'GB ' R ' ^RR' RR' R" • • • > V 1Z J 
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where the subscript ^ stands for the contribution to 0(p n ) chiral constant. For our further 
discussion we will explicitly need 

4b = x (u ^ +x+) - (13) 

The leading order of the odd intrinsic parity sector of Cq b coincides with the Wess-Zumino- 
Witten Lagrangian [35] Cyyzw' ^ or * ne ex P nc it form of even parity part £q B and complete 

Cqb see IS EE], ( see also P3)- 

The most general interaction Lagrangian which is relevant for the saturation of the 
0(p 4 ) LECs [23j is linear in resonance fields, namely 

4 4) = c d (Su^} + c m {S X+ ) + id m (P X -} + i^(P)( X -} + 



1 /2 

and all the couplings are of the order 0(N C ). This is true also for the last term of the first 
line with two traces which is enhanced due to the rj exchange (see Appendix [S] esp. (|112p ). 
This term with d m Q (depending solely on the singlet component of P) has not yet been studied 
in the phenomenology as it always contributes to the saturation of LECs together with the 
large Nc enhanced r\ exchange. The complete operator basis of the 0(p e ) even intrinsic 
parity of R^T has been constructed in [21]. 

Integrating out the resonance fields at the tree level we reconstruct the Lagrangian C x pt 
of ChPT, schematically 



exp J d A xC x PT^j = J f-Rexp J d A xCn x T^j ■ 



(15) 



Effectively up to the order 0(p e ) the integration over R is equivalent to the insertion of the 
solution R( 2 } of the lowest order equation of motion (i.e. those derived from £^ Rkin + £^ B ) 
for resonance field R into the Lagrangian Cr x t- Because the resonance fields R couples 
to the 0(p 2 ) building blocks in and the resonance masses are counted as 0(p°), we are 
consistent with the chiral counting (jllj) . Finally we get 

^ = 4) +4) +4) + ... 

with explicit separate contribution from Goldstone bosons part of the R/\T Lagrangian Cqb 
and the leading Nq contribution of the resonances 

4 2) = 4?b (i6) 

4*> 2 ) = + (17) 



where particularly 



r (4) _ (A*) , ~(4)\ I 

S:,_R - \ L 'RR,kin~*~ i -R J \R-+R(?) 

r (6) _ f r (6) r (6) r (6) r (6) 

^X,R ~ \'~RR,km~ t R RR' RR'R" I \r~^R( 2 )' 
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The structure of Lagrangians C GB and C R are identical to C x , just the couplings are 
different. Then for generic chiral coupling constants k x of £ x pt we may write 

k x = k GB + k XtR , (18) 

where k Xt R corresponds to the resonance contribution. The usual hypothesis of resonance 
saturation assumes k G B to be very small and the resonance contribution k X) R is expected to 
be dominant. 

The above resonance saturation strategy and the construction of all relevant operators 
were studied already in the past. In [23] was found the basis for all relevant resonance 
operators contributing to 0{p i ) and their contribution to LECs while in [24] the authors 
presented the extension to 0(p e ) in even- intrinsic-parity sector. 

In this paper, we complete this effort presenting the construction of basis and resonance 
saturation at 0(p 6 ) in the odd-intrinsic parity sector. 



3 Lagrangian of RxT in odd-intrinsic parity sector 

Before starting the construction of resonance monomials let us summarize the structure of 
the pure Goldstone-boson part of the odd-intrinsic sector. The leading order starts at 0(p 4 ) 
and the parameters are set entirely by the chiral anomaly. The Lagrangian is given by |35| 
(see also [8]; note we have the same convention for the Levi-Civita symbol, i.e. eoi23 = 1) : 

CwZW = 4^ e ^{/ ^^4^) -i(W^a()(U,l,r) , (19) 

with 

^ nvafi = L^LyLaRp + jL^RjyLaRp + \L^ v L a Rp + \Rp U L a Rp — \a pL v R a Lp + OpR va Lp 
— cr^UvRaLp + a fl L u L aj 3 + a fl L l/a L j 3 — i<jpL v L a Lp + \op_LyOa_Lfi — ia v a a L p — (L -<->• R) , 

where we have defined 

L fl = u l^v) , L^u = ud ll l v u\ R^ = uW^u , R^ u = u d^r v v) , Op = {v) , dpu] 

and (L <-> R) stands also for a •<-)■ a* interchange. The power £ indicates a change of u 
to v£ = exp(i£<j)/(Fy/2)). Concerning the 0(p 6 ) part we will stick to the form introduced 
in [8]. Let us only note that we will drop the index r and the explicit dependence on the 
renormalization scale \x from the corresponding LECs C^, but one should have in mind that 
any studied in this text is a renormalized LEC with the scale set to some reasonable 
value (~ M p , M^i) to make good sense of the following study. 

For the construction of the operator basis in the odd intrinsic parity sector of RxT we 
use the same tools as in |24j . where the reader can find further details. First we construct 
all possible operators built from chiral building blocks and resonance fields that are invariant 
under all symmetries. Then in order to find the independent basis we use 

1. Partial integration 

2. Equation of motion 

V%4(x— ^<X->) (20) 
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3. Bianchi identities 



V/JVp + V„r p „ + VpT^ = for = - y +flu (21) 



4. Shouten identity [36] 

9crp^a/3fiu Qaa^fi^vp Qofi^-pvpa 9ap,^-vpafi ~\~ Qav^pafip = (22) 

5. Identity 

V V = " 2 K * r /-] " V "/-^ (23) 
All relevant operators in odd parity sector can be written in the form 

Of = eT^df^ , (24) 

with the basis for individual monomials Of a p, with X = V , A, P, S, VV, AA, SA, SV, 
VA, PA, PV, VVP, VAS, AAP; so Lagrangian becomes: 



X i 



The basis of the operators Of^ ua/3 is summarized in Tables HE We have included there only 
the operators relevant in the leading order in the 1/Nc expansion i.e. operators with one 
flavour trace and those with two traces that are enhanced by rj exchange (see Appendix lAl 
for details). This represents main result of our work. 

As is shown in |24| 133]. we can further modify the resonance Lagrangian (|25p . The reason 
is that the resonance fields play merely the role of the integration variables in the path 
integral (|15p and can be therefore freely redefined without changing the physical content of 
the theory. As a consequence we can choose appropriate field redefinition in order to eliminate 
some subset {Of}^ x ,i)eM °f 0(p 6 ) operators from Cr^^ an d shift their influence effectively 
to the 0(p 6 ) terms including the remaining operators {Of}(x,i)dM an d also to the higher 
chiral order terms C^^° , symbolically 

4 6 X T dd) = £^c^ £ ^ f+4>v dd ) (26) 

(X,i) (X,i)fM 

The possible new terms £^r° dd ^ of the order 0(p 8 ) and higher generated by such a redefini- 
tion can be neglected because they do not contribute to the 0(p & ) LECs when the resonance 
fields are integrated out. Note however, that after such a truncation we get new Lagrangian 



4 6 X T dd) = £ (27) 

{X,i)(£M 

which is not equivalent with the previous one on the resonance level. On the other hand, the 

LECs obtained form £^,y dd ^ coincide with those derived form C^^. . 

The stars in the Tables HE indicate those operators which can be eliminated by the 
resonance fields redefinitions discussed above and means therefore a redundance of a given 
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i 


b v R 


i 


O v R 


1 


\{V^ u {h a(T u u u 13 - u /3 u a h aa )) 


11 


\ v \J + J J — j /ypcr 


2 


\(yt lv {u a h a ' T u li - uPh aa u a ) 


12 


\ v \J + ' J /iipcr 




\(V^ u (v n l3 h aa — h^v^v )) 

1\V \LLq-Uj IV It tli l*i (J J I 


1 3 


1 \ v J+ MA. — / 


4 




14 




5 




15 


i(V^[/5,X+]) 


6 




16 




7 


Uy^luaf^vf* ~ U P f« a uJ)) 


17 


(^{v ff /r/}) 


8 


i{V> xv {f^ T u a u ts - uPu a f°: a )) 


18 


(X^<i/V)(X-> 


9 


(V^{ X -,u a uP}) 






10 


(V» u u a x-u p ) 







Table 1: Monomials with one vector resonance field. 



i 


QA 

i fivafi 


i 


i [iva.fi 


1 


{A^[u a u^ U(T u a }) 


10 


\{A^u a ){V^ X -) 


2 


{A^ u [u a u a u^u a }) 


11 


i(A^{f^, X -}) 


3 


(A^{V a h^,u a }) 


12 


i(^{V«x-X}> 


4 


i(A^[ft P ,u°u a }) 


13 




5 




14 


i(A^[f^, X+ }) 


6 


\{A^{fl u u^u a - u a uPfl a ) 


15 


(A^{V a £%u a }) 


7 


[{A» V (u a f^uP -vPfSfUa) 


16 


(A^{V a f^,uV}) 


8 


{A» u {fT ,hP a }) 






9 


i(A^f^)(x-) 







Table 2: Monomials with one axial-vector resonance field. 



monomial as far as its contribution to the resonance saturation is concerned. The details are 
shown in Appendix [Bj Note, however, that this redundance concerns only the saturation and 
not the direct calculation of the correlators and amplitudes with resonances in the initial and 
final states. We will return to this point later on. 

In the following section we will demonstrate the use of the resonance basis for two classes 
of examples. The resonance saturation will be studied in Section [5j 

4 Applications 

In this section we illustrate applications of the Lagrangian usm S t wo examples. We 

study two three-point correlators, namely (VVP) and (VAS), and use both OPE constraints 
as well as phenomenological inputs to fix the relevant coupling constants. In the first case we 
also discuss some phenomenological applications in more detail. 



10 



i 










/H5 














1 


(P{f 


4 


{Pu^u v u a u^) 


1 


{S[f a ^u^}) 


2 


i(Pu a ,/fV) 


5 




2 




3 


i(p{/r,« a ^}) 


6 









Table 3: Monomials with scalar or pseudo-scalar resonance field. 



i 


Operator Of^ Q/3 , i? = y,^ 


Operator O™ a0 ,R = P,S 


1* 


i{R^R a ^){X-) 




2* 


i({R^,R a P}X-) 




3 


({V CT i?^,i? QCT }^) 




4 


({V^i?^,i? aCT }n a ) 





Table 4: Monomials with two resonance fields of the same kind. 





Operator 


i 


Operator 6f%, a0 


Operator 


1* 


i(L4^,S]/f) 


1* 


i([V» u ,S}fT) 




2* 


u a uP]) 


2 


i([V' iV ,V a S\vf i ) 





Table 5: Monomials with two resonance fields of different kinds. 



i 


Operator 6™ a0 


i 


Operator O[^ a0 


i 


Operator Of^ va0 


1* 




1* 


({A^,P}f^} 


1* 


i({V^,P}u a uP) 


2* 


i(V» u (A aa u aU P - u^u a A acT )) 


2 


{{A^,V a P}uP) 


2* 




3* 


i(V^(A ar7 uPu a - u a u^A aa ) 






3* 


{{V^,P}ff) 


4* 


\{V^{u a A aa u^ - u^A aa u a )) 










5* 


({V^,A a r>}f^)g pa 










6* 


i{[V^,A^] x +) 











Table 6: Monomials with two resonance fields of different kinds. 



i 


Operator 6?™ 


1* 




2* 


i([V» v ,A a P]S) 


3* 


{A^A^P) 



Table 7: Monomials with three resonance fields. 
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4.1 VVP Green function revisited 

The standard definition of this correlator is 

K°(P^) = f d 4 xd 4 ye i ^-y(0\T[V-(x)V, b (y)P c (0)\0) 

with the vector current and the pseudoscalar density defined by 

A a A a 
V^(x) = q( x ) lfl —q(x) , P*{x) = q(x)i l5 — q{x) , 



(28) 



(29) 



(our convention is 75 = Z7°7 1 7 2 7 3 ). This correlator was already studied in the past, see e.g. 
[TT] . [12] . [28] . Here we provide a complete result based on our £r x T) i-e. also with two- and 
three-resonance vertices that were not considered in [28]. Using Ward identities and Lorentz 
and parity invariance one can define 



H(p) J? = ^^ntp 2 ,? 2 ; r 2 ) . 
The OPE constraints dictate for high values of all independent momenta 



n((A P ) 2 ,(Ag) 2 ;(Ar) 2 ) 



BqF 1 p 2 + q z + H 
2A 4 



+ ° 1 * 



5 F 2 1 „ / J 



whereas in the case when only two operators are close to each other one gets 

n((Ap) 2 ,(g-A P )V) 

n((Ap)V;((/ + Ap) 2 ) 



A 2 p 
1 1 



2.,2+°[ A3 



(30) 
(31) 

(32) 
(33) 



In the following we will use only two-momentum OPE. The reason is that for general cor- 
relator, not all the high energy constraints can be simultaneously satisfied using only finite 
number of resonances in the effective Lagrangian. This statement has been proved in [T7] for 
the case of the {PPS} three-point function. For the (VVP) this problem has been partially 
studied in [12] (see also [37]). 

By means of an explicit calculation based on (the relevant Feynman graphs are 

depicted in Fig(T]) we get 




® 





Figure 1: Feynman graphs contributed to VVP Green function. Double lines stand for 
resonances and dash lines for GB (double lines together with dash lines is the sum of both 
possible contributions). The crossing is implicitly assumed. 
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N c 4F v K^ v p 2 \Qy/2d m F V K PV 32d m K P 



167r 2 r 2 r 2 (p 2 - M 2 )(q 2 - M 2 ) - M 2 )(r 2 - M 2 ) r 2 — Mp 



qj p2^yyp 9 p2 

+ ^2 _ M 2)U _ A/f 2 \ L 8K 2 ^3 J 



(p2-M2)(g2-M2)(r2-M2) (p 2 — My){q 2 — My) 

[ P 2 (kY 6 + 2k\ 2 ) - q 2 (i4 6 - 2k\ 7 + 2k\ 2 ) - r 2 (^\ A + k& + 2k&)] 



2V2F V 



r 2 (p 2 - M 2 ) 

+ (p^q). (34) 



From OPE (J3TJ we get then the following constraints for the couplings 

V _ N C v , v N C v N C p _ n 

KU ~ 256V2« 2 F V ' Kl6 + 2Kl2 ""32^F v ' * 17 " " My/2**F V ' ^ " ' 

vv _ F 2 + 16y/2rf m F v «f y jV c M 2 yy yy_ F 2 

2 ~ 32F 2 5l2n 2 F 2 ' ^ 2 " 3 " 8^" ' [6b) 
By employing these constraints one getH 

1 F 2 C p 2 i 2 , 2\ y c4 1fW p2^VVP 

1 T-rRyT/ 2 2 2^ _ f ^ g 1 4tt 2 F 2 100^^ 

S W >« > J 2 (p 2 - M 2 )(g 2 - M 2 )r 2 (r 2 - M 2 )(p 2 - M 2 )(q 2 - M 2 ) 

16V2d m F v 4 v [(p 2 + q 2 )M 2 - 2r 2 M 2 ] 

r 2 {r 2 - M 2 ){q 2 - M 2 ){p 2 - M 2 ) ' 1 j 

This should be equivalent with the LMD+P ansatz introduced in [11] so that two independent 
constants b and c introduced there are directly connected with phenomenological couplings 
K PV and k vvp . Considering just vector resonance interactions, k vvp = k pv = (or equiv- 
alently taking the limit Mp — > oo in (|34p ). we can reconstruct the LMD ansatz [12] 



1 p2 (J2 , n 2 , „2\ N C M y 

J_ n LMD f J 2. 2x _ £1 . CP + ? + r J ~ 4^ ( o 7) 

Bo *>*> T >- 2 (p2 _M 2 )(g 2 -M 2 )r 2 - (3?) 

The result in ChPT up to 0(p 6 ) at the leading order in 1/Nc expansion includes two LECs 
from the 0(p e ) anomalous sector 

i-n— (p 2 ,q 2 ;r 2 ) = -^+32^ - + g 2 ) . (38) 

Comparing this with a low energy expansion of the RxT result (|34p we give the following 
lowest-resonance contribution to Cj and (cf. also Section [5]) 

w = F 2 (84 y -4 V ) , d m F 2 K VVP _ V2d m F V K PV 2d rn K P _ F v {2k\ 2 + 8< 4 + 
7 8M* + 2M 2 M4 MpMy + M 2 4^2M 2 

22 " V2M 2 2M4 ' 



2 Note that these constraints imply automatically also the fulfilment of (|33[). However, the requirement (|32[1 
cannot be satisfied until k^ v = 0, which is in contradiction with another high-energy constraint for related 
pion transition form factor; see next subsection (cf. also |12]). 
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Using the OPE constraints (|35p we obtain 
CV = —rrr + 



F 2 d m F v {-2^14 v M 2 + F v k vvp ) 
7 64M V + 2M 2 M y 

rW _ , N C _ 2V2d m F v 4 v 

22 16M y + 64ir 2 M 2 My ' [ ' 

4.1.1 Formfactors 

Let us define fully off-shell formfactors for V*^*^* vertex, where V can represents either pion 
(or any other Goldstone boson) or pseudoscalar resonance via 

T Vll (p\q 2 ;r 2 ) = -^( r 2 -m 2 v )U(p 2 ,q 2 ;r 2 ), (41) 

where Z factor interpolates between pseudoscalar source and V. Let us discuss in detail the 
7r°77 formfactor. We have Z n o = 3/2BF and using the OPE constraints (1350 we can define 
(note we are working in the chiral limit) 

^ R o*>V;r 2 ) = ^r 2 n R * T ( P V;r 2 ) , (42) 

where II RxT (p 2 , q 2 ; r 2 ) was introduced in (|36p . For on-shell pion the k vvp drops out (note 
that this is not connected with the chiral limit simplification) and we get a simple result 

P („2 , „2Wi , QQ./o d m F v pyx Ng M v 
^TT^ ' q ' Uj " 3 (p 2 - M 2 )(q 2 - M 2 ) ■ [6) 

.PV 



Dropping Kg we can again reconstruct the LMD ansatz 



F w P + 



2 i n 1 _ N C M 



^o 77 (P ,9,0) k ^ =q - ^o 77 (p , C? , 0) - T (p2 _ M 2 )(g2 _ M 2) • ( 44 ) 



Using Brodsky-Lepage (B-L) behaviour for large momentum [38] [39 

1 

Q 

one can arrive to the following constraint 



B-L cond.: lim .7^0^(0, -Q 2 ; ml) ~ — — ^ , (45) 



7?2 

B-L cond.: k PV = = . (46) 

32V2d m F v 

Before discussing the possible violation of the Brodsky-Lepage condition let us study the 
influence of the constraint (|46|)on the original VVP Green function. The n R * T correlator 
in (|36p will now depend only on one constant k vvp and we get 

1 ttRxTV 2 2. 2, _F 2 (P 2 + q 2 +r 2 )-^l F 2 [{p 2 + q 2 )M 2 - 2r 2 M? 



v\ 



— Il R x T (p 2 q 2 -r 2 ) =— KF ^- j A^W + 

B ^ ' q ' ' 2 (p2 - M 2 ){q 2 - M 2 )r 2 2 r 2 (r 2 - M 2 ){q 2 - M 2 ){p 2 - M 2 ) 
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1 fW F 2 k vvp 

{r 2 - M 2 ){p 2 - M 2 ){q 2 - M 2 )' 1 j 

The violation of the OPE (|32p is here manifest. The remaining constant k will drop out 
for an on-shell pion in the formfactor and one gets: 

B-L cond.: T R ? T (p 2 , q 2 ; 0) = ^ ■ (48) 

^TT^ ' y 12vr 2 F(p 2 -M 2 )(q 2 - M 2 ) y 1 

For completeness let us also provide the ChPT result. From (|38[) using (|40p we have 

J2 , „2 , ^2 772 



x-ChPT/ 2 2. m _ 



l + ^(l-i^(l + 32v^ 



k pv- 



(49) 



and 

B-L cond, ^V,^ ) = -^(l + ^). (50) 

For a reader's convenience let us also summarize previous results based on the VMD and 
LMD+V ansatzal : 



rVMD, 2 2. m _ N C My My 

^ ^ ' 9 ' ^ 12^ (P 2 - My) (g 2 — My) ' 



rLMD+V/J „2. m = ^ p 2 g 2 (p 2 + g 2 ) + /^(p 2 + g 2 ) 2 + h 2 p 2 q 2 + /l 5 (p 2 + q 2 ) + h 7 

[P,q,) 3 (p 2 — M Vi )(p 2 — M 2 2 )(q 2 — M Vi )(q 2 — My 2 ) ' ^ 

with (valid in the chiral limit) 

7 " 4vr 2 F 2 ■ 
We have therefore the following relation (compare with (144 p ) 

7$%tf, t\ 0)| B L = ^(P 2 , q 2 ; 0) . (53) 

Now let us turn back to Brodsky-Lepage condition. We have seen it has important conse- 
quences on the actual form of the 7r° — 7 — 7 formfactor within RxT. However, recent BABAR 
measurement [40J showed phenomenological disagreement with this condition. There are also 
theoretical arguments |17] which showed that high-energy constraints cannot be all satisfied 
for a given formfactor within the ansatz with only finite number of resonance poles. (For 
a recent study on Brodsky-Lepage revision see [5T]; see also [12] and references therein.) 
We will thus relax the Brodsky-Lepage condition by allowing a small deviation from (|46p 
parameterized with 5bl 

K^, -iv^bv <1+^BL, • (54) 

Its actual value can be set by fitting the BABAR and CLEO data. In this fit and also in the 
following phenomenological applications we set 

M v = m p K, 0.775 GeV, M P = rrv ~ 1-3 GeV, F = F w ss 92.22 MeV (55) 



3 The LMD+V ansatz adds one extra vector multiplet in comparison with the LMD one. This corresponds 
to MHA for which all the OPE and B-L constraints can be satisfied simultaneously. 
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Figure 2: CLEO (blue points) and BABAR (green squares) data with fitted function 
•T^o* (0) — Q 2 ;0) defined in ([4~3"|) using the modified Brodsky-Lepage condition in (f54"|) . The 
full line is for 5bl = —0.055 and (blue) dotted line stands for standard B-L (i.e. 5bl = 0). 
Dot-dashed (red) line shows fitted function A(Q 2 / (10 GeV 2 ))^ with A = 0.182 ± 0.002 GeV 
and /3 = 0.25 ± 0.02 as obtained by BABAR collaboration [40J. The asymptotic 2F is repre- 
sented by the horizontal dash line. 



and also (for details see [43] ) 

F v = F p = 146.3 ± 1.2 MeV, d m ps 26 MeV . (56) 
The new BABAR data indicates important negative shift in <5bl with the result 

5 BL = -0.055 ± 0.025 . (57) 
Our fit together with CLEO and BABAR data is depicted in Fig. [2j 



4.1.2 Decay p 



7T7 



In this subsection we illustrate a particular phenomenological application of the above results, 
namely a prediction for p — > ivy decay. For this process we can use a connection with the off- 
shell 7T77 formfactor introduced in the previous subsection. First, let us define the amplitude 
A for the process p + {p) — > ir + (p)'y(k) (we will use only the charged decay process to avoid 
the discussion on to — p mixing for the neutral one): 



pol. P 



(58) 



16 



from which we have already factorized out the Levi-Civita and momentum dependence. Sim- 
ilarly one can define the amplitude for tt°(p) — > 7(^)7(2) 

T n o^ = ^~J2 \^^ ua0 k a lpe;(k)el(l)\ 2 — . (59) 

pol 



The connection with -^77 formfactor is obtained via 

Ap^ry = I™ a (q 2 - My)T „- 77 (0, q 2 ; 0) (60) 

I by My q2^My 

and quite simply 

An— » 77 e 2 Jvn(0,0;0). (61) 
Putting these two definitions together we can extract the ratio and corresponding parameter 

2eF v 



Ap^Tfy 



A 



1 + x . (62) 



>77 

Using the experimental value Tp^ n ^ = 68 ± 7keV this parameter was obtained to be equal to 
x = 0.022 ± 0.051 in [11] based on the 1992 edition of the particle data book (same number 
was also used later e.g. in [12"]). Updating this prediction with a new experimental input we 
can get flip in the sign 

exp: x = -0.003 ± 0.054 . (63) 

The change is mainly due to a new value of Fy (study e.g. in [43]) and a new precise 
measurement of ir° lifetime by PrimEx group [H] (see also [30] ) . 

Within our formalism, the parameter x defined above is proportional to the deviation from 
the simple VMD ansatz (|5ip . or in other words from the exact Brodsky-Lepage condition (cf. 
(j53j) ). Using (j43j) and (|54l) we get in terms of 5bl 



4vr 2 F 2 

X = MPc feL • (64) 

The results of the previous subsection allows us to make rather precise determination of this 
value 

R X T: x = -0.010 ± 0.005 , (65) 
which using ([62]) and experimental input for Y^o^^ leads to the following prediction: 

RxT: Y p ^ wl = 67.0 ± 2.3 keV . (66) 



4.1.3 Decays of ?r( 1300) 

In the previous section we have obtained a prediction for the p — > ttj decay width. However, 
it was based on the ratio of two decay widths (cf. (|62p ) and experimental input of one of 
them. We could predict also the absolute value for p —> irj directly from (|58p and (|60p 
without the necessity to use the experimental value of tt° — > 77 (in fact we will discuss a little 
the latter process in the very next subsection) but one should remember that we have been 
making several simplifications, namely: we are working in large Nc, using only lowest-lying 
resonances and we are in the chiral limit. All together within this approximation we cannot 
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expect the accuracy of the result being better than 30%-40%. On the other hand one can 
expect that some of the systematic uncertainties will cancel out in the ratios similar to one 
studied in the previous part. 

The same strategy can be repeated for 7r(1300) decays. In fact we can work in exact 
correspondence; the two decays would be now: 7r(1300) — > /try and 7r(1300) — > 77. The only 
problem now is that none of these two processes have been seen so far. The most recent limit 
on 7r(1300) — > 77 by Belle collaboration [45] 

IV^ Tr <72eV (67) 

sets at least rough limit in our studies. Using the definition ([41 p the main object here is 



Tp rt & ' q ' m P) ~ 3 Fv (f - M^)tf -Ml) ■ (68) 

The amplitude for 7r(1300) — > 77 is given by 

An'-^-y-y = e 2 J r p 77 (0, 0; m 2 /) (69) 
and similarly for 7r(1300) —> (see also (|60|) ). Then 

RxT 2 8V2 gygK^Mg - F v k vvp 

,RxT 4,/2 y^TMg - FyK™ 

3M V M v • 1 ' 

Both these amplitudes depend on one so far undetermined constant k vvp . Provided we have 
experimental values of both branching ratios we could verify the consistency of our model. 
In the present situation we can visualize how one decay mode depends on the second one, 
and this was done in Fig. [3J One can see that we have two solutions for k vvp as we have 
quadratic equation for decay width as a function of k vvp and none of these two solutions 
can be ruled out. Note that the full width for 7r(1300) is assumed to be between 200 and 600 
MeV (see [IS]), so both processes are extremely suppressed for any of these two solutions. 

The experimental bound on IV can be used to get estimate of k . In order to fulfill 
the limit (f67[) we expect the numerator in ([68|) to be suppressed. This expected suppression 
leads in analogy with (|54l) to the following ansatz 

F 2 M 2 

™' = -l6^ (1 + ii) ' (72) 

with parameter 5a which should be reasonably small. In terms of this parameter we get the 
decay width in a compact form 

and thus the extreme phenomenological suppression of 7r(1300) — > 77 can be understood 
within our formalism to be due to the small factor (5bl — ^a) 2 - The experimental limit 
together with (|5T|) set the allowed range for the parameter 5a 

- 0-27 < 5 A < 0.16, (74) 

which is good enough to set the value of kvvp to 

kvvp ~ (-0.57 ±0.13) GeV. (75) 



18 




Figure 3: The connection of decay width for 7r(1300) — > 77 and 7r(1300) — > pj (note that 
we have two possible solutions). The dashed line represents the Belle's limit (|67|) on T^i^m 
(grey area is thus excluded by this experiment). 



4.1.4 Decay ir° — > 77 and r\ 



77 



As we have stated, the absolute decay widths are accessible via our approach only with the 
limited precision. For instance for the 7r° — > 77 amplitude we have obtained only very simple 
prediction (|6ip . It turns out, however, that it agrees very well with the experimental determi- 
nation. On the other hand, similar determination for r/ — > 77 would be a phenomenological 
disaster. 

In order to go beyond the leading order we can use the chiral corrections calculated using 
ChPT. The most recent study of 7r° — > 77 amplitude went up to NNLO [30]. The motivation 
for going beyond NLO lays in the fact that there are no chiral logarithms at NLO \47\ 148] . At 
NNLO these logarithms though non-zero are relatively small so the play very important 
role. We can therefore use existing calculations within ChPT with our estimate (|40p of 
C| . Here our approximation, namely the chiral limit, does not make any difference as by 
construction LECs [Cff in our case) do not depend on light quark masses. With previous 
phenomenological determination of the couplings k? v and Kyyp we obtain 



°"' " MM* + 2 W P ^ - * ( °- 35 ± ° M} X 10 ~ 



" 3 GeV" 2 . (76) 



The second and last unknown LEC at NLO for ir° — > 77 and r\ — > 77 is . Anticipating 
the result of the next section and using the OPE constraints (|35|) we get 



c w = N c N c _ F v kY 3 F v k\ v d m0 F 2 dmoFyKwp 

8 768M 2 vr 2 512tt 2 M2 ^2M 2 2M^ + md m M 2 p M 2 QMpMy ' 1 ' 

where we have also dropped the term proportional to 5bl because it is not numerically 
relevant. Unfortunately at this moment similarly as for already mentioned d m Q we cannot 
make an estimate for n( 3 and kY (all these couplings are dominated by the rf exchange, cf. 
Appendix [A]). We may however again connect two-gamma decay widths of 7r° and rj. We 
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may for example set the unknown from the experimental value of T{r} — > 77). This was 
done for NLO rj — )■ 77 expression in [30J. There is ongoing project which should enlarge this 
calculation to the NNLO within ChPT (for preliminary results in the chiral limit calculation 
see [49]). We thus rather postpone as a future project the final determination of 7r° — > 77 
based on the experimental value IV^™. Let us only mention, that if we assume that the 
NNLO corrections for rj are indeed small as for ir° — > 77, the value in (|76p has roughly the 
influence at 0.5% level for r„.o — > 77 (with the opposite sign). A new study of isospin breaking 
effects in ]50J indicates another shift of the similar size (however now with a positive sign) 
and thus at this moment we do not expect quantitative change of the prediction made in [30 . 

4.1.5 g-2 

Probably the main motivation for studying the VVP correlator is hidden in the determination 
of the muon g — 2 factor. It is beyond the scope of this paper to discuss this problem in great 
detail (for details see e.g. [51] ). Let us only mention that the main source of theoretical error 
for its standard model prediction comes from hadronic contributions, more precisely from 
the hadronic light-by-light scattering which cannot be related to any available data. The 
hadronic four-point correlator WW is further simplified into three classes of contributions: 
a) tt^ and loops b) 7r°, 77, 77' exchanges and finally c) the rest, which is modelled usually 
via constituent quark loops. It is clear that this separation is not without ambiguity and 
different approaches can differently calculate contribution especially between a) and c) or b) 
and c). Our contribution based on the VVP correlator study belongs to the group b). To 
avoid inconsistency we will work in close relation with similar work done for LMD or VMD 
ansatze [52] . Using the fully off-shell (i.e. including also the 7r° off-shellness) tt° — 7 — 7 
formfactor (I42D we arrive to 

a LbyW> = ( Q5 8 ± L2 ) x 1Q -11 _ ( 7g ) 

In the error only the uncertainties of our model were included. The systematic is mainly 
influenced by the above mentioned ambiguity of how one defines and splits the pion-pole and 
regular part from the (WW). We have put the cutoff energy at 10 GeV. For a better 
comparison let us present in Table [8] predictions for the studied contribution to the muon 
g — 2 for the different models summarized in (I52p . We have recalculated there the light-by- 



model 


a]l hyhy ' x 10 11 


VMD 


57.2 


LMD 


73.7 


LMD+V "on-shell" 


58.2 


LMD+V "off-shell" 


72 + 12 


this work 


65.8 ±1.2 



Table 8: Contribution of ir exchange to the muon g — 2 factor for different models. 

light contributions based on VMD and LMD ansatze. We have also reevaluated the case 
of LMD+V ansatz or more precisely its on-shell simplification as defined in (|52p . Three 
unknown constants are set similarly as done in [12], i.e. hip = and /15 is based on the 
p — > 7T7 phenomenology /15 = 6.99 (obtained for the updated value in ([63]) ). The full LMD+V 
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"off-shell" ansatz has 7 parameters (for details see [37] ). One relation can be obtained from 
the chiral anomaly and others can be: Brodsky-Lepage behaviour, higher-twist corrections 
in the OPE and one-large momentum OPE, together with data (CLEO for this turn) we are 
still left with two undetermined parameters. Their variations in reasonable range set the final 
error for the corresponding LMD+V value in Tab. [5J Let us note that also the possibility of 
B-L violation together with new fit of two parameters (hi and /15) was studied in [37] with 
no influence on the central value of g — 2 contribution. Too many parameters is not the only 
problem connected with the LMD+V ansatz. Status of p(1450) as a first radial excitation 
of p(770) is doubted by the possible existence of lighter p(1250) [53] • Its presence is also 
supported by the study within AdS/QCD approaches [M]. On top of that the inclusion of 
the complete set of all excitations in all channels (i.e. inclusion of it") can change again the 
studied ansatz similarly as we have encountered for the first excitation (see (|44|) and (|53p). 
Let us also note quite astonishing coincidence of our result with the most recent study 



based on AdS/QCD conjecture [55] af = 65.4(2.5) x 10 
4.2 VAS Green function 

The (VVP) Green's function studied in the previous section represents without any doubts 
the most important example of the odd intrinsic sector of QCD. However, it is not the only 
quantity one can analyze using our complete lowest-lying resonance model. As the second 
example we have chosen (VAS), which has not yet been studied (to our knowledge) in the 
literature. It also enables to demonstrate the use of the "second half" our the odd intrinsic 
resonance Lagrangian, i.e. those with 1 ++ and ++ states. 
Defining (beware of the same symbol as for (VVP)) 

<f(p,<?) = J d 4 xd 4 ye i P- x+ ^y(0\T[V-(x)A b u (y)S c (0)\0) , (79) 

with (cf. also ([U])) 

A a A a 
Al(x) = q(x)~i^ b —q(x), S*{x) = q(x) — q(x) . 

Similarly as for VVP one can write 

n(p, q)f v c = f ab %^p a qV,q 2 ;r 2 ) , (80) 
where r = — (p + q). In resonance region, we have 

1 tj/ 2 2 2, 8V2F y «-2< 5 ) 16y/2F A K& , 32c m 4 , 16V2F ACm Kf A 
-U(p ,q -r ) = -= ~o ^72- + To ^72 + 
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B ^ ' y ' ' p 2 -M 2 q 2 ~M\ r 2 - M 2 (q 2 - M 2 A )(r 2 - M 2 ) 

8V2F v c m (2 K f v + 4 V ) \QF a F v kX a 16F A F v c m K VAS 

(p 2 — My)(r 2 — Mj) (q 2 - M\)(p 2 - M 2 ) (q 2 - M\)(p 2 - M 2 )(r 2 - M 2 ) ' [ ' 

At high energies one can obtain the following OPE relation 

n((Ap) 2 , (Xqf- (Xrf) = B f d 2p2 ~t~/ +o(^] (82) 



2A 4 p 2 q 2 r 2 \X e 
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and again we will not consider here one- momentum OPE limits. The high-energy constraint 
leads to 

F 2 



K 2 



c 14 — ' 
TP fO^SV , „SV\ 



2k V 
ZK 15 I 



,VA 



32F A F V ' 



2F A ref 



SA 



16^/2c 



If we use these relations, we have finally only one free parameter: k vas ; the result is 



-DO 



M 2 +M\ + Ml) + 32F A F v c m K VAS 



(83) 



(84) 



2{q 2 - M 2 A ){p 2 - M 2 )(r 2 - Mj) 

From the theoretical point of view we are thus in a better position than we were for (VVP). 
After imposing OPE we are left with one free parameter whereas in the case of (VVP) we 
had two (cf. ([36]) ). We can thus simply connect all processes schematically represented as 



(V:p,u,K*, -,,...) ~ {A: ai ,h,K x ,GB,W...) ~ (S : a,K,a ,f ,H . . .) 



(85) 



via a single parameter. The problem is that they are very rare and have not yet been measured, 
on top of that the status of some of the particle content is controversial by itself (especially 
if talking about a scalar sector). The parameter k vas can be, however, set using other (not 
that rare) processes it enters. One way is to check in next section to which of 23 parameters 
it contributes and use directly the system of LECs. This can be done already here within the 
calculation of VAS. At low energies, up to 0(p 6 ) one has 



1-U(p 2 ,q 2 ,r 2 

-DO 



-2,2CYi 



(86) 



Comparing with the low energy expansion of the full RxT result (|8ip we get 



F A K 



V2M 2 



it + M* 



2< 5 ) 



A 

F A C m K^ 



2V2M 2 



+ 



c m n§ , F a F v kY a , c 



Ml 



+ 



2M\Ml 



+ 



2V2MIM 2 



y/2M\M 2 s 



F A FyC m K VAS 

2M\M 2 S M 2 



(87) 



Using the OPE (jHSJ) we obtain 

64 



M 2 M 2 



+ 



M\M 2 



M\M 2 



VAS 



+ 



F A F v c m K 
2M\M 2 M 2 



The knowledge of Cyy leads directly to the value of k vas and thus to the rare processes 
schematically specified above. The current attempts for a Cy[ estimation were summarized 
in Table 1 of [56j with rather inconsistent values obtained both from the phenomenology 
([57], |58j ) and by a model-dependent determination |56| . The most precise value seems to 
be obtained from K+ -> data [59]: Cf x = (0.68 ± 0.21) x 10" 3 GeV" 2 [57]. Using the 
values set in ([55|) and (f56|) . together with 



m 



"0 



984.7 MeV, 



42 MeV 



and the Weinberg sum rules (to get values of M A and F A ) we arrive at 



VAS 



0.61 ± 0.40 GeV. 



(89) 
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4.3 Short note on the field redefinition 

The previous two examples were calculated using the full resonance Lagrangian £^,^, dd \ 
Here we would like to address a question what would happened if one would repeat the same 

calculation but instead use the reduced resonance Lagrangian 4|° dd) . This L agrangian is 
established in Appendix [B] and can be obtained from the full Lagrangian f)25 j) by means 
dropping the operators marked with a star in the Tables [T][7] i.e. by means of omitting 20 
parameters: k rr , ^f A , K>f v , k-Y A , k pa , k pv , k rrr and using the bar over the rest of nf (see 
Section TB.4p . This can be motivated by its equivalent contribution to the saturation of LECs. 
This exercise was already performed in [23] for (VAP) with an interesting finding, that after 
imposing the OPE condition the both results are the same. In our case the conclusion is, 
however, different. Using the reduced resonance Lagrangian we would not be able to simply 
fulfill the OPE constraints by imposing some conditions on . In the first case, the OPE 
for (VVP) requires an additional relation, namely My = . In the second case, (VAS), 
the OPE cannot be satisfied at all. 

Thus we have to conclude that the equivalence of both calculation in the even sector for 
(VAP) was just a coincidence and it is not a general feature. 



5 Resonance contributions to the LECs of the anomalous sec- 
tor 

We have seen in the previous two applications the explicit examples of the calculation with 
the resonance fields. A match between this result in a region of small momenta (i.e. p <C Mr) 
at one side and the ChPT result at other side enables to extract the dependence of LECs on 
resonances. In this way we have obtained within VVP calculation and (see (I39p ) 
and from VAS it was possible to extract Cy[ (|H7)1. The dependence of all others on 
the parameters of the resonance model can be obtained by systematic integration-out of all 
resonances. So obtained Lagrangian can be expand over the canonical basis of NLO odd- 
intrinsic Lagrangian established for example in [8]. In this way we have saturated 21 of 23 
constants and only and C^g stayed intact as they are subleading in large Nc- The rf was 
explicitly considered (see Appendix lAl and [60] ) and it contributes in Cg , and C^. It is 
always the first term in these LECs and we put it in the boldface font to stress its large Nc 
dominance over the rest. Generally we have the following expansion in large Nc for all , 
schematically 

Cf = ai N c + kNc + 0(N C ) , (90) 

where a% 7^ for i = 6, 8, 10 and 6j = for i = 3, 18. 

The field redefinition similarly as done in |24] was performed and details are summarized 
in Appendix|Bj All 20 parameters denoted by stars in Tab {TH 7]can be dropped in the following 
and for all others a bar should be added (bar parameters nf are defined in the last section of 
Appendix [B]). We prefer, however, to use the original parametrization as it represents direct 
connection with the resonance phenomenology and is thus simpler to use. 

The explicit form of the resonance saturation generated by the resonance Lagrangian (|25p 

is: 

r w _ d m 4 2V2d m G V K\ v _ V2d m G V Kl v V2G V 4 _ 2^2G v k% V2G v kY 

1 M 2 p M 2 p My MpMy My My My 
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My 2My MpMy ' 

w = FaK^ C m Kf C m F A Kl A y/^CmGyKf 1 C m Fy ft| V | F V k\ _ y/lGyK^ 

2 ~ V2M\ M| ^2M\M 2 M 2 M 2 2^2M 2 M 2 2^2M 2 M 2 

F a G v kJ a K WAS c m F A G v 
M\M 2 M\M 2 M 2 ' 

C™ = 0, 

yW _ F a k% d m K% V2d m G v n PY F a F v kJ a V2F a k£ F A nf 2 F A nf 5 



4 2V2M 2 M 2 MpMy 2M\M 2 M\ ' 2^2M\ 2V2M 2 
F\k aa d m F A n PA d m F v n PY 2c d 4 ^/2c d F A K SA V2c d F v 4 v 

AM\ 2\[2M 2 A Mp V2M 2 M 2 M 2 M\M 2 M 2 M 2 

| c d F V Kf | F v k\ V2F v 4 F v k% F v 4 F v 4 | V2G v kY 4 
V2M 2 M 2 V2M 2 M 2 y/2M 2 V2M 2 sf2M 2 M 2 
F A F v kJ a F a F v kJ a _ 2F v G v k yy K VAS c d F A F v _ K YYP d m F v G v 
M\M 2 2M\M 2 M v M\M 2 M 2 M 2 M y 

w _ F A nf 2 _ d m n% d m F A n PA _ d m F V K Py F v kX q _ G v k\ 7 _ F v G v k yy 

5 ~ V2M 2 A M 2 J2M\M 2 V2M 2 M 2 ^J2M 2 ^J2M 2 M y ' 

w ^ N c , F A 4 _ F a F v k ya _ F v G v k yy _ 2V2F a k£ F AK f F AK f 5 

6 576M2vr 2 3^2M\ 3M\M 2 3M y 3M\ ^2M\ 3^2M\ 

F\4 A 4c d /cf 2^2c d F A K\ A 2^2c d F v 4 Y V2c d F v 4 Y V2F v k\ 
6M A 3M§ 3M\M 2 3M 2 M 2 3M 2 M 2 3M 2 

2y/2F V K% V2F v k v V2F v k\ y/2G V K^ 3 G v k\ 7 F v k\ & 
+ 3M 2 + 3M 2 + 3M 2 + M 2 3^2M 2 ~ ^2M 2 

2F a F v k ya F a F v k ya 2F v G v k yy 2 K YAS c d F A F v V2F A d m0 K, PA 2d m0 K% 



3M\M 2 3M\M 2 My 3M\M 2 M 2 3M\M 2 3M 2 

dmQK2_ d m0 F V K PY \/2"d m0 Fy/tf v 2^d m0 GyK PY _ 2d nl0 F V GyK YYP 

3M 2 3V2M 2 M 2 3M 2 M 2 3M 2 M 2 3M 2 M V 

^ w _ 2d m K,£ y/2d m F V K PV F v k{ 2 V2F v kY 4 F v k\ & 



QW __ _ _ 

7 M 2 ~ ' M 2 M 2 2V2M 2 ' M 2 ' A^/2M 2 



F 2 K VV ??2,,VV K VVP^ F 2 

My 8My 2M 2 My 



c w N C F v k\ 2 _ F v kY 3 F v kY s F 2 kY V F 2 kY 

8 768Mg7r 2 6V2M 2 ^2M 2 \2^2M 2 2M y 2AM y 



V2d m0 F v n py d m0 F 2 K YYP 2d m0 n p 



3M 2 M 2 6M 2 M V 3M 2 ' 

nW _ F A 4 F 2 k aa F a k£ V2F a k^ F a i4i F a4 5 

l^o — 



9 4V2M 2 8M A 2V2M 2 M 2 " \[2M 2 A AV2M 2 
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F A K 2^ j_ 2d m Kf \f2d m F A K^ A d m F A 4 A k aap d m F^ 
~M| MjM 2 V2M\M 2 p 2M\M 2 ' 

cW = N C ^4 ^14 A , _ Fa4 F a k? F a k? 5 

10 768Mg7r 2 \2^2M\ 2AM\ 6V2M\ \f2M\ 2^2M\ \2yf2M\ 

F\^i A F 2 K AAP d m0 V2F A d m0 K PA F A d m0 n PA 2d m0 K? 
2M\ 6M\M 2 3M]M p 2 Z^2M\M% 3M 2 ' 

w _ F A nf A c m n% c m F A nf A c m F V K,f v c m Fy«;| v F v k% F v kY 5 
— — /= T H TTo 7= — n o H f= — n — T H ;= — n — T + 



11 \/2M| M| v/2M^M| \/2MjM y 2^/2MfM y 2^/2M y 
F a F v kJ a K YAS c m F A F v 



2M^M y 2.U 2 ,.U;.U l - 
w _ ^/2Gy< _ V2gwg _ y/2gyg Gy4 V 

12 ~ Ml M 2 Ml + M y 

w ^ F a k^ F A F v nJ A FyGyn^ 2V2F AK j F a k* _ F 2 k aa Ac d n s 2 

13 " V2M] Miikf2 Ml + M| V2Mj 2Mj + M| 

2^/2c d F A K^ A 2V2c d F V Kf V | ^C^Fy/^ | ^Fy/tj 7 2^/2FyK^ \/2FyK 



V 



+ 



M\M 2 S M'iM'l M 2 s Ml Ml Ml Ml 

V2Fy4_ V2G v kY 2 G v k\ & V2G v kY 7 2F a F v kX a F a F v k^ a 
Ml Ml v2.U f ' + Ml + M\Ml + M|M y 

2K VAS C d F4Fy 



w = F a k a _ F v 4 F a F v kJ a V2F AK j F AK f 5 _ F 2 k aa 2c d 4 
14 " 2\~ni\ \~2Ml 2M\Ml M\ 2^2M\ \M\ M 2 

V2c d F AK SA V2c d F V Kf y c d F V Kf F v k\ ^2F v k\ F v k\ F v k v 



+ : r o 'n + ^ ' „ + 
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M\M 2 M 2 Ml V2M 2 Ml ^ M y \/2M£ ^ 

F a F v kY a F a F v kJ a K VAS C d F A Fy 

M\Ml 2M'\\\'l M\M 2 Ml ' 

w = F a k a F a F v kJ a _ F v Gv4 v , V2F A Kj F a k a _ F 2 k aa 2c d 4 

15 " 2^2M\ 2M\Ml Ml M\ 2^2M\ \M\ + M 2 

V2c d F AK f A y/2c d F V Kf y c d F v 4 Y , F v k\ ^2F v k\ F v k\ F v k, 



V 



M\M 2 M 2 Ml ^2M 2 Ml ^M 2 M y ^2 My ^2M y 



gy^ F A FyK x A F A F V 4 A K VAS C d F A Fy 

V2Ml M\Ml + 2M\Ml + M|M|M2 ' 

w V2F a k,2 F A Gv4 A _ y/2Gv4 , vggyg _ yggyg F A G v /cjf A 

16 M] + M]M2 M y M y Ml + M|M2 

w = F A 4 _Cd4__ V2c d F A 4 A V^CdFyKfV _ V2c d G v 4 v F a F v k( a 

17 \2.\l 2 M 2 M\M 2 M 2 Ml M'z-M'1 M\M% 
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F a G v kJ a Fa4 V2F a k£ F AK f 5 _ F 2 k aa 2c d 4 c d F A n s 2 A 



M\M 2 V2M A M\ 2^2M\ AM% M 2 s/2M 2 A M 2 

c d F v 4 w y/2Fv4 V2G v k\ F v k% F v 4 F a F v ^ a F a F v kJ a 

2 V / 2M 2 M 2 Ml Ml ^2Ml V2Ml 2M\Ml 2M 2 A Ml 

FvGvk£_ G 2 v kJ v F v Gv4 w K YAS c d F A F v K VAS c d F A G v 

AMI + 2M V 4M V + M\M 2 s Ml M\M 2 s Ml ' 

°18 — U ! 

w= y/2F AK £ F a G v kJ a V2F V 4 yggy^ F y < 6 G v 4 6 V2G v k\ 7 

19 M\ M\Ml Ml Ml J2MI \/2My. Ml 
F A F V 4 A FyGyK™ F 2 kJ v F v G v kY v 

+ M\Ml + My 2My + My 

w = V2F A 4 _ F 2 k aa _ F 2 kJ v FyGyK™ _ F a k£ _ F a k£ F a k? 5 2qk| 

20 Ml 2M\ AMI 2 My \f2M\ V2M 2 A ^2M\ M 2 
V2c d F AK f A V^CdFynfV c d F v 4 v V2F V 4 s/2Fyi$_ \/2FyKg 

M\M 2 M 2 Ml j2M 2 M 2 Ml Ml M? 



A lV±g ^S ±V± V V±±Vl s ±V±y lV±y lV±y IVAy 

V n..~V jp.jp.. ~VA jr.rn.~VA z7.z?..,,VA jP,.n..~VV 

3 



Me _ GWg FaF v 4 a F A F V 4 A F a F v 4 a _ FyGyK 
2V2MI V2Ml M\M 2 M 2 M 2 M 2 M 2 My 
K VAS c d F A F v 



+ 



M\M 2 Ml 



w _ V2F A 4 F 2 k aa F a F v kJ a F 2 k vy F v Gv4 y , Fa4 F AK f 5 



2 



21 " M\ 2M\ 2M\Ml AMI 2M V V^M 2 y/2M_ 

2c d 4 V2c d F A 4 A _ V2c d F v 4 Y _ c d F v 4 w y/2F V K% F v k% j F v 4 



F v 4 F y < 6 gy^ GyKyj F A F V 4 A F A F V 4 A F A F V 4 A 
V2M 2 2V2MI V2MI y/2Ml M\Ml 2M\Ml 2M\Ml 

3 



F v Gv4 y n YAS c d F A F v 



2M V M\M 2 Ml 



F,rK V F 2 K VV 

qW _ r v K n _ r y K 3 
22 V2Ml 2M V ' 

jp.^A jp2 ^AA 

n W _ r A^\% r A K 3 / 01 \ 

623 ""7iMi"^Mr' ( } 

Apart from already mentioned relations = and C|g = we have found out one further 
relation free from nf 

P V r<W _ tp (nW r<W\ , nff 



, ( , v = F V (C& -C&) + G V C™ . (92) 

The transformation established in Appendix [B] was employed as an independent check of the 
previous relations. 
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6 Summary 



In this paper we have studied the odd-intrinsic sector of the low-energy QCD. We have 
constructed the most general resonance Lagrangian that describes the interactions of the 
Goldstone bosons and the lowest-lying vector-, axial-, scalar-, pseudoscalar-resonance multi- 
plets. We were working in the large Nq approximation and considered only those terms that 
contributes to 0(p 6 ) anomalous Lagrangian (i.e. to the first non-trivial order). This was the 
main aim of our work. We then demonstrated the use of this Lagrangian for three different 
applications. The first two represent calculations of two three-point Green functions (VVP) 
and (VAS). The third application was the complete integration out of the resonance fields 
and establishing the so-called saturation of LECs by resonance fields. 

The first application VVP is the most important example of the odd-intrinsic sector, 
both from the theoretical and phenomenological point of view. We have discussed different 
aspects of this Green functions. First, after calculating this three-point correlator within 
our model and imposing a certain high-energy constraint we ended up with the result which 
depends only on two parameters. These were further set using new BABAR data on ^77 
off-shell formfactor and Belle collaboration's limit on tt' — > 77 decay. After setting these two 
parameters we can make further predictions. The outcome of our analysis is for example a 
very precise determination of the decay width of a process p — > ttj: T p ^ n ^ = 67(2.3) keV. 
We have also studied a relative dependence of the rare decays %' — > 77 and tt' — > pj. Based 
on the experimental upper limit of the former one can set the lowest limit of the latter. 
Prediction of our model is 30 keV > T^i^.^ > 4 keV (based on Belle's IV-^ < 72 eV). 
Next, we have also evaluate the value of LEC together with short discussion on 7r° and 
r\ two photon decays. Last but not least a very precise determination of the off-shell ir°- 
pole contribution to the muon g — 2 factor was provided. Our final determination of this 
factor is = 65.8(1.2) x 1CT 11 . The RxT approach has thus reduced the error of the 
similar determination based on lowest-meson saturation ansatz by factor of ten and is in 
exact agreement with the most recent determination based on AdS/QCD assumptions. Let 
us note that the present theoretical error for the complete anomalous magnetic moment of 
the muon is around 50 x 10~ n and the experimental error around 60 x 10~ n [61] (with the 
well- know discrepancy above 3d). A new proposed experiment at Fermilab E989 [62 J plans 
to go down with the precision to the preliminary value 16 x 10~ n and thus the reduction of 
the error in the theoretical light-by-light calculation is more than desirable. 

If VVP represents very important and rich phenomenological example, the three-point 
correlator (VAS) is connected with very rare processes and represents so far never studied 
example of the odd-sector. We have established its OPE behaviour which enabled us to reduce 
the dependence of the VAS Green function to one parameter. This opens the possibility of 
a future study of these rare but interesting processes. 

In the last section we have studied the resonance saturation at low energies. We have 
integrated out the resonance fields to establish the dependence of LECs of odd-sector C t - on 
our parameters. As we are limited by large Nc we cannot make prediction for C3 and C}g 
but we have set all other 21 LECs. We have found one relation among , C^, and 
C22 free from our parameters. 
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A The large Nc counting 
A.l General considerations 

Let us start with the Ul(Nf) x Ur(Nf) invariant Lagrangian for the nonet of the GB and 
resonances without using the equations of motion and the Cayley-Hamilton identities. Then 
the large Nq behaviour of the couplings accompanying individual operators in the effective 
Lagrangian with octet GB (after rj' has been integrated out) can be understood as follows. 
Let us write in the same way as in 



u = e^ T °/ F ^ 2 u, (93) 

where T° = y/l/N F l and 



e i<f> a T a /FV2 



is the SUl(Nf) x SUr(Nf) basic building block, and therefore 



^° = fv ^Mdetn). (94) 

Let us also remind |24j . that (ft and (ft a do not mix under the nonlinear ly realized Ul(Nf) x 
Ur(Nf) symmetry. For the construction of the Ul{Nf) x Ur(Nf) effective Lagrangian, we 
have the usual building blocks constructed from u and the usual external sources lu,,ru,x an d 
X + (now also with singlet components) e.g. 



F 



1 / 2 .o. 



etc. at our disposal. In the above formulae, the covariant (in fact invariant) derivative of 0° 
is defined as 

iV° = <V>° "2a°F, (96) 
however, it does not represent an independent building block because of the identity 
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The above set of building block have to be further enlarged including also the external sources 
9 for the winding number density 



W = 4^ trcG » vG ^ ' (98) 

with covariant derivative 

= d^O + 2a° . 

We have to include also the following invariant combination 

dP 

X = 9 + ^-. (99) 
F 

Let us remind the large Nq counting for the generating functional of the connected Green 
function of quark bilinear s and winding number densities 

Z[l,r,x,X + ,0} = N^Z [9/Nc}+N c Z l [l,r,x,X + ,e/Nc} + ... , (100) 

where the ellipses stay for the subleading terms in the 1/Nc expansion. This implies the 
usual Nc counting of the physical amplitudes with g glueballs and m mesons 

A g , m = O(N 1 c +Sm0 - g ~^). (101) 

This counting should be reflected within the construction of the effective chiral Lagrangian 
of R X T. 

According to the (|101|) . the explicit resonance fields have to be counted as 0(N C ). 
As far as the GB are concerned, within the tilded building blocs, each member of the pseu- 
doscalar nonet is automatically accompanied with (minus)one power of the decay constant 

1/2 

F = 0{N C ' ), which ensures the right counting of the vertices with GB, provided the corre- 
sponding fields are counted as 0(N C ). The only subtlety is connected with the field cjP . 

The origin of the field (jp in the individual terms of the Lagrangian is twofold. It can either 
come from the tilded building blocks Y = u«, h^ u , x± (and from their covariant derivatives 
D^Y ; note that it completely decouples from and ftf?) or from the X— dependence of the 
Lagrangian. Each operator O constructed form the tilded building blocks only (and therefore 
including at least one flavour trace, the only exception is O = 1) is in general accompanied 
by a potential Vq{X) which is a function of the variable X only, 

C = Y,V d (X)d. (102) 

6 

While (jP originating from the tilded operators is counted as 0{Nq) as the other GB, however, 
the same field coming from the power expansion of the potentials counts as 0(1/Nq) within 
the large Nc expansion. Therefore, expanding the general operator O and the corresponding 
potential Vq(X) in powers of (jP and its derivatives (and taking into account that F = 

1 12 

0{N C )) we have the following natural rule for the order O(N^) of the resulting coupling 
constant at a term of this expansion with T flavour traces, R resonance fields and no fields 

2-T-^R-^n <n<2-T-^R-^n . (103) 
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The lower or higher bounds are saturated in the case when all </>°'s come exclusively either 
from Vq(X) or from O. 

Suppose that we had used the LO GB equations of motion prior to the expansion in 
powers of cjp. This allows to eliminate the terms with derivatives, namely [23] 



Such a transformation of the original tilded operator do not create any extra trace in contrast 
to the octet case. Because the singlet mass Ml = 0(l/N c ), the cfP dependence of the 
resulting operator brings about a factor of the order 0(N C , 3 ^ 2 ) (the same, as if (fP came from 
the potential) and the above bounds on n remain therefore valid. On the other hand, the 
further simplification using the Cayley-Hamilton identity can destroy them, provided we use 
it in order to eliminate terms with less traces in favour of the terms with more traces. 

The next step is to integrate out (fP treating the mass M5 2 as O(p ). This can be done using 
its equation of motion, derived from the corresponding part of the LO Lagrangian expanded 
in powers of (fP 

Cf = l -D<tP ■ Dtp - \Ml {<tPf - i-|=( x _)/ + d (P)^ + . . . , (105) 

where do term comes from the expansion of the potential and is therefore of the order 0(Nq 1 ). 
The solution for i^ reads in the leading order of the p expansion^ 



= 0{N 3 C !' ! ) + O(JVg), (106) 

where we have depicted the orders of both terms, (jp^ should then be inserted into the 
original Lagrangian expanded in powers of cfP . As a result, taking ()103[) into account, the 
orders of the multiple trace operators within the SUl(Nf) x SUr(Nf) operator basis are 
enhanced. Namely, we have the following bound for the corresponding couplings 

13 11 

2 - T - -Ro - -n {P) < n < 2 - T - -Rq + n (x _> - -n (P) , (107) 

where Tq and Rq are the numbers of the traces and resonance fields before elimination of (fP 
and n/ x _\ and are the numbers of the new factors (x~) an d (P) (which appear after cjp 
is integrated out) respectively. More conveniently this can be expressed in terms of the actual 
number of traces T = To + + n^ x _^ and resonances R = Ro + as 

2 — T - -R + n (x _ } <n<2-T-^R + n {P) + 2n (x _ } . (108) 

The loophole of this formula is, that for its application one has to trace back which of the 
factors (P) and (x—) originate in the <jP dependence of the tilded Lagrangian. The extreme 
cases are either none or all of them, which gives a much raw estimate 

2-T-^R<n<2-T-^R + N {P) + 2iV (x _) , (109) 

where now N( x _} and are the total numbers of (P) and (x-) traces in the operator, the 
lower bound corresponds now to the usual trace and resonance counting. 



'Here we have took into account, that the resonance fields should be counted as 0(p 2 
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A. 2 Explicit examples 



Let us illustrate the above statements by means of an explicit examples. For instance, the 

— 1/2 

coupling at the term (Sx—)\X—)i a ^ first sight of the order 0(N C ) might be of the order 

1/2 3/2 

0(N C ) or even 0(N C ), because it can originate either from the term 



i(Sx-)W {3x _ } (X) = i{S( X - + ...)) (™{Sk-> X + ■ ■ ■) 



, (HO) 



which has the constant w}g x \ = 0(N C 1 ^ 2 ), (this corresponds to the lower bound (1108P ) or 
from the term 



(Sx + )W (Sxt) (X) = (S (x + - j-^/x- +...)> (w? s , +) + ■ ■ ■ 



where «;? s , = 0(N^/ 2 ); (this corresponds to the upper bound (|lU8p ). 

Similarly the coupling d m o at the operator i(P)(x-) (see (fT^j) ). naively of the order 
0(N G 1 ^ 2 ) can be enhanced by the (f>° exchange. Indeed, inserting ()106|) to the term do{P)(/> 
of the Lagrangian (I105p . we get the following contribution to d m o 



dp FVN^ _ nrAT i/2, 
Ml 2V2 



dmo = ^k-^ = O(N^) , (112) 



where we have taken into account that do = 0(N C , 1 ). 

Let us give also some examples of the odd intrinsic parity terms with resonances, which 
similarly to the previous example lead to Nc enhanced multiple trace terms when cfp is 
integrated out. Some terms with one resonance are for example 

Cr = w(^15^^])^ra(I) + e^(F7f)^(I) + ^(#7+Vra(^) , (113) 
where 

W m (X) = Y. w m xk ( 114 ) 



and where 
These generate the operators 



ffi=0, w V=0(N c 1/2 ), for . = 1,2,3. (115) 



6 V n = is^(V^ff)(x-), 
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6$ = i(A^ff)( X -) (116) 

1 /2 

with the couplings of the order 0{N ( j ) (i.e of the same order as analogous single trace 
operators and therefore included in our basis) and 

8^{V^[u a ^]){P), 

{A ^ff)(P) (117) 

with the couplings of the order 0{N^, 1 ) suppressed with respect to the single trace operators. 
The two-resonance example is 

Crr = e^{V^V aP )W RR1 (X) + e^{A^A a P)W RRl (X) , (118) 

where 

W RRl (X) = Y. w Rm xk with ^ = °> ™ i £ ) m = 0(Nc 1 ), for z = 1,2. 

k 

It gives rise to the operators 

6t A = ie^(A^A^)( X -} (119) 

with the couplings of the order O(Nq) (the same order as the analogous single trace operators 
and therefore included in our basis) and 0(N C ) operators 



e^(A^A^)(P), (120) 

which are suppressed with respect to the single trace ones. 

As the last step, we integrate out the resonance fields in order to get the resonance 
contribution to the odd parity sector LECs of the resulting \PT Lagrangian. It can be done 
using the 0(p 2 ) EOM for the resonance fields and inserting their solution back to the 
RxT Lagrangian. The general form reads 

R(2) = _L jg) (121) 

Mj, R y ' 

where J R 2 ^ = 0{p 2 ) comes from the LO resonance Lagrangian flHD- Because J [ R 2) = 0(N l J 2 ), 
the order of the contribution of the individual terms of the RxT Lagrangian (with <jp inte- 

1 /2 

grated out) can be obtained counting the resonance fields as 0(N C ). This gives finally the 
following simple bound on the order of the contribution of the operator with T traces, total 
N/p\ factors (P) and total -^v"( x _) factors (x~) originating in the to the LECs 

2 - T <n<2-T + N {P) + 2iV (x _ } . (122) 

The lower bound represents the usual trace counting. Note however, that the upper bound 
have to be taken with some caution, because it can be saturated only in the case when all (P) 
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and (x~) traces appear as a consequence of the cfP dependence and that this ^dependence 
comes solely from the tilded operators and not from the potentials. For a given operator these 
two conditions need not to be satisfied simultaneously. 

The fact that the Nc of order some operators can be enhanced could further complicate 
the usual way of the saturation of the ChPT LECs. Namely, in the process of integrating out 
the resonances, it is assumed, that loops can give only NLO contribution suppressed by the 
factor 1/Nq for each loop. This counting could be apparently complicated by the enhanced 
operators. Let us illustrate this point assuming the contribution of the following term of the 
odd RxT Lagrangian £ 

I = ... + W 2 AP (X)e, ual3 ({A^,V a P}^) + ... 

= ...-2 W i P e^{A^P)^J-^ + ... (123) 
with = O(Ng). This gives rise to the following enhanced Nc term 

- ir F W * P M* W<^V"P)^( X _) = 0(N C ) . (124) 

Apparently, this term contributes to the 0(p s ) LECs, when the resonances are integrated out 
at the tree level. However, the bubble with two such vertices gives a contribution to the 0(p 6 ) 
operator d a (x-)d a (x-) of the enhanced order 0(N C ). The same is true also for analogous 
operators from the even sector, e.g. 

VnXHiD^P}^) = 2vf p {PD,S)^-^f- + ... (125) 

with vf p = O(Ng) which leads to the enhanced operator 

iv^ P ^< P VWx-> = 0{N C ) (126) 



counted as 0(p 8 ) in the tree level saturation process. The bubble with two vertices 

i"' P W (p8 " s>8 " (x - > <127) 

leads to the expression 

ir {ir F v " P ^i) S Ad V(x-w)s v (x-w) 

d d k jHx - y) d d p PtxPv 



e ih-(x-y) ~ r r/irv /-^ 

(2Tr) d (2vr) d (p2 - Mj + iO) ((p - k) 2 -Mp + iO) V ; 
+0( P 8 ) 

and (after addition of appropriate counterterm) results in the following 0(N C ) contribution 
to the coupling Cd a (x-)d a (\-) associated with 0(p 6 ) operator d a (x~)d a (x-) 

r PS-loop _(^_\ V Z_ 2 > V Z_ ' (-1<>Q\ 

u d a (x-)d-(x-) ~ 64vr 2 VM V Ml -Ml ' { ' 
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Though the above loop contribution are enhanced by the factor Nq with respect to the 
naive trace counting, it does not mean, that loop counting fails. The reason is that the LO 
contribution to CQi x _)d a (x-) that comes from the tree level and originates in the kinetic term 
of the field / 

\W - ~\ {Wfm) '*<*->*•<*-> = ( Nl) (130) 

so that the loop contribution is suppressed by l/Nc as usual. 

Let us finally comment briefly on one point, which also might lead to confusion. In [24J, 
the following operators are abandon using the large Nc arguments, namely 

i(Pu^)( X -), 
i{SP){ X -), 

i(V^ X -)(P)- (131) 
These can be, however, derived from the operators (before doing any transformations) 

i(Pu^)V(X) , 
i(SP)V(X), 

i(V M V"x + )F(X) , (132) 

by means of integrating out the field 0°, which appears from the potential for the first two 
operators (and saturates therefore the lower bound of (|1Q8|) ) and from the building block x+ 
for the last one (and corresponds therefore to the upper bound of (|108p ). According to our 
rules the operators are of the order 0(N^/ 2 ), O(Nq) and 0(N^/ 2 ) respectively (as the similar 
operators without additional trace) and all of them contribute therefore at the 0{Nq) order 
of the LECs of the effective chiral Lagrangian staying at the operators 

(x-v M )(x-}, 
(x+x-)(x-) , 

(V^V"x-)(X->- (133) 
However, these operators can be derived analogously as above from 

(X+X+) , 

(V^x+) , (134) 

by the process which saturates the upper bound of (|108p and results in the order O(N^). The 
abandoned operators lead therefore to the NLO contribution to the corresponding LECs. 

B Field redefinition 

As we have discussed in detail in Section [3j by means of appropriate field redefinition we can 
effectively eliminate subset of the 0(p 6 ) operators from the Lagrangian C^'^^ and shift their 
influence on the ChPT LECs into the effective coefficients nf which stay at the remaining 



34 



operators of the chiral order 0(p G ) and higher. As a consequence, the 0(p 6 ) LECs resulting 
from the process of integrating out the resonance fields from the Lagrangian Cr x t depend 
only on these effective couplings k£ which are particular linear combinations of the original 
resonance couplings nf . In order to identify these relevant combinations and the redundant 
operators, we can proceed in several steps. 

B.l Elimination of 0\J , O aa , O VVP and O aap 
With the field redefinitions 

v -+ v - w («r (x-)v^ + i^ v {x-, v#} + l r vvp {p, ^}) , 

we get for the 0(p 4 ) part of the Lagrangian 

/•(*) , -(4) ~(4) -(4) 

l -'RR,Hn '"R ~^ *~RR,kin "r ^fl 

- /-< 4 ) _j_ 

~~ 1 -RR,kin R 

-k\ v oX v - 4 v oX v - k vvp o vvp - 4 A o AA - 4 A o AA - k aap o aap 



v^M^V 1 13 2 ° 14+ 2 K ° 3 
v^M 2 



v^M 2 



+ k^o£ + X -k aap O p A ^j + 0(p 8 ) . 

At the same time, the same redefinition applied to £^,£ dd ^ generates only the additional 
terms of the order 0(p 8 ) and higher, which can be neglected as described above. We can 
thus eliminate the operators 0\\ ', , O VVP , O aa and Q AAP and include their influence on 

the 0(p 6 ) LECs effectively into the constants kY 3 , kY 4 , k pv , kY 8 , Kg, k pv and k a , k a 1: k pa . 
B.2 Elimination of Of A and O VAS 

In the same way we can eliminate also the mixed bilinear terms using the field redefinition 

My \ 

+i4 A (A afS u a up - upu a A afS ) + i4 A (u /3 A Q/ V - uM°%) 
+i^(u CT u !y V m/3 - V a ^u u u a ) + i/^ A (u CT F Q/ V - u v V a ^u c ) 
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We get then 



- K por A - 4 A o 2 yA - 4 A o 3 yA - 4 A or A - 4 A o 5 yA - 4 A o\ A - n VAS o VAS 

and the operators Oj A and O^^ 5 are thus eliminated. The only relevant additional effect of 



the redefinition comes from transformation of C 



(4) 
R 



Fy 
F A 

iGy 
2^2 



(V/D 



2^ 



-^ a o 4 a +4 a (o 6 a -^4 a ) 



+ 4 A (O s A - ^O a ) + k\ a O a - ^ a O a + k 



.VA/nA 



.VASqSA 



2^2 
iGy 



2^2M\ 

K 



,va^v 



va^v 



va^v 



Gy 



2V2M? 



V 



-2k\ a O a + 4 a (O a - Of) 

.VA {0 A_ A ) + K VA (p A_ A ) 

+ 4 A (0 A - A ) + 24 a O a + 2k vas 0| A 



2^/2 

+ 4 A {o A -o A ) vx\-\oi - 



Here we have used that 

({A^,fl p }f^)g^e, uaa = 
and other similar consequences of the Shouten identity. 

B.3 Elimination of Of A , Of v , Cf y and Of A 

Finally we can further eliminate another terms by the redefinitions 



S 
P 



A 



[IV 



V, 



{XV 



S + ^ W (i«f A [/f > A n + >4 A [u a u0, A^}) , 

V - ^ W (}4 A [S, ff] + k S 2 A [S, u a u?] + «f A {P, 

V - ^ w (i«f ^} + i4 v u<*p u e + 4^{p, /f }) . 



We get then 



-1m|(PP) - ^Mf <SS) + iMi(^V) + \M V {V^V, V ) 
l -M 2 P {PP) - l -Ml{SS) + \m\(A^A, v ) + ±M V (V^V»„) 
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-K PA PA - K PV PV - K PV PV - K PV PV - K SA SA - K a/i O a/i - K aV O a 
1 ^1 ™2 2 3 3 1 1 2 2 1 1 

therefore the operators Of A , Of v and Of ^ are eliminated. We get additional contributions 

c m {S X+ ) -> c m (5x + ) + ^(-4 A Of4 + 4^3-4 y Or 5 ) 



B.4 The effective couplings nf 

Putting the result of previous subsections together we get the parameters of the repa- 

rameterized and truncated Lagrangian , which is relevant for the saturation of ChPT 

LECs, as a functions of the parameters nf . As we have discussed above, the LECs have to 
depend on the couplings nf of the original Lagrangian only through their particular 

combinations . We have proved this by means of direct calculation as a nontrivial check 
of the formulae (I9ip . 

< = k\ 

^2 — ^2 
K 3 = K 3 



K V K V °d f SV | F A VAS \ F A _y A 



,.\ - K V _ F A ,.\ A 



2\/2M 



2 n 3 



v v f a va 

y F A y A 

s 2v^Mi 2 



v 2Gv4 v _ ( pv _ 2G v k vvp 

Kg + V2M V 2Ml P 1 2^2M2 
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V u m PV 

Kl ° ~ 2M 2 " 2 

v I Fa va 
Kl1 + 2~V2M 2 ^ 




2V2M^ K4 



Kg 



A , d m Q I PA Fak aap \ Fak aa 

Kg + 



K A 



N F M 2 \ 1 2^2M 2 J V2M 2 A 



Fak aa d m ( PA Fak aap 



^2M 2 2M 2 \ 1 2V2M 2 



K A 
K 12 
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"13 



"14 



"15 



"16 



K 2 



K sv 



K 13 + 



K 14 + 



Kl5 
K 16 



«f + 



G v 



,VA 



+ 



L v 



Fy 



.VA 



2V2M 2 



2M| 



2M| 



,SA 



,SA 



V2My' 

F v 



.VAS 



,VAS 



2\/2M£ 



Fa 



,SA 



2V2M 2 A 



+ 



K 2 



F A 



2V2M 2 A 



F A 



.PA 



Gv 

Fy 
2^2 M? 



.VAS 



+ 



Gy 



.sv 



V2M? 



+ 



F A 



.VAS 



K 



VAS 



F A K 



V 
AAP S 



2\[2M\ 



2\[2M\ 



4 + 



Fv 



.PV 



2V2M 2 1 



l v 



K 3 



F v 



„PV 



2G v k 



VVP 



2y/2M v 



2y/2M v 



G\ 



V2My 



2GyK 



VVP 



2V2M 2 



p 

K 5 



4 y 

K2 



Fv 



2V2M V 



K 3 V 



FyK 



VVP 



2V2M V 



+ 



V 

Fv 



2\[2M 2 A 



,SV 



2V2M 2 



+ 



F A 



,VAS 



2V2M 2 A 



Gy 



.PV 



FyK 



VVP 



V2M 2 



K 2 V 



2V2M 2 



References 

[1] S. Weinberg, Physica A 96 (1979) 327. 

[2] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

[3] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 



[4] J. Bijnens, Prog. Part. Nucl. Phys. 58 (2007) 521 [arXiv:hep- ph /0604043 1 . 

[5] J. Bijnens, G. Colangelo and G. Ecker, JHEP 9902 (1999) 020 [arXiv:hep-ph/9902437] . 

[6] J. Bijnens, G. Colangelo and G. Ecker, Annals Phys. 280 (2000) 100 
[arXiv:hep-ph/9907333] . 



39 



[7] T. Ebertshauser, H. W . Fearing and S. Scherer, Phys. Rev. D 65 (2002) 054033 
[arXiv:hep-ph/01 10261] . 

[8] J. Bijnens, L. Girlanda and P. Talavera, Eur. Phys. J. C 23 (2002) 539 
[arXiv:hep-ph/0110400l . 



[9] J. F. Donoghue and E. Golowich, Phys. Rev. D 49 (1994) 1513 [arXiv:hep- ph/9307262|. 

[10] M. Davier, L. Girlanda, A. Hocker and J. Stern, Phys. Rev. D 58 (1998) 096014 
[arXiv:hep-ph/9802447] . 

[11] B. Moussallam, Phys. Rev. D 51 (1995) 4939 [arXiv:hep-ph/9407402] . 

[12] M. Knecht and A. Nyffeler, Eur. Phys. J. C 21 (2001) 659 [arXiv:hep-ph/0106034| . 

[13] P. Masjuan and S. Peris, JHEP 0705 (2007) 040 [arXiv:0704.i247l [hep-ph]]. 

[14] M. Golterman and S. Peris, Phys. Rev. D 74 (2006) 096002 [arXiv:hep-ph/0607152] . 

[15] G. 't Hooft, Nucl. Phys. B72 (1974) 461. 



[16] S. Peris, B. Phily and E. de Rafael, Phys. Rev. Lett. 86 (2001) 14 [arXiv:hep-ph/0007338] . 



[17] J. Bijne ns, E. Gami z, E. Lipartia and J. Prades, JHEP 0304 (2003) 055 
[arXiv:hep-ph/0304222] . 

[18] I. Rosell, P. Ruiz-Femenia and J. Portoles, JHEP 0512 (2005) 020 
[arXiv:hep-ph/051004T]. 



[19] I. Rosell, J. J. Sanz-Cillero and A. Pich, JHEP 0701 (2007) 039 [arXiv:hep-ph/0610290 |. 



[20] A. Pich, I. Rosell and J. J. Sanz-Cillero, JHEP 0807 (2008) 014 [ arXiv: 0803 .15671 [hep- 
ph]]. 

[21] I. Rosell, P. Ruiz-Femenia and J. J. Sanz-Cillero, Phys. Rev. D 79 (2009) 076009 
[arXiv:0903.2440l [hep-ph]] . 

[22] A. Pich, I. Rosell and J. J. Sanz-Cillero, JHEP 1102 (2011) 109 [ arXiv:1011.577T1 [hep- 
ph]]. 

[23] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[24] V. Cirigliano, G. Ecker, M. Eidemuller, R. Kaiser, A. Pich and J. Portoles, Nucl. Phys. 
B 753 (2006) 139 [arXiv:hep-ph/0603205] . 

[25] J. Portoles, AIP Conf. Proc. 1322 (2010) 178 [arXiv: 1010. 33601 [hep-ph]]. 

[26] E. Pallante and R. Petronzio, Nucl. Phys. B 396 (1993) 205. 

[27] J. Prades, Z. Phys. C 63 (1994) 491 [Erratum-ibid. C 11 (1999) 571] 
[arXiv:hep-ph/9302"246] . 

[28] P. D. Ruiz-Femenia, A. Pich and J. Portoles, JHEP 0307 (2003) 003 
[arXiv:hep-ph/0306i57] . 

40 



B. Ananthanarayan and B. Moussallam, JHEP 0205 (2002) 052 [arXiv:hep-ph/0205232] . 



K. Kampf and B. Moussallam, Phys. Rev. D 79 (2009) 076005 |arXiv:0901. 46881 [hep-ph]] . 

S. R. Coleman, J. Wess and B. Zumino, Phys. Rev. 177 (1969) 2239. 

C. G. . Callan, S. R. Coleman, J. Wess and B. Zumino, Phys. Rev. 177 (1969) 2247. 

G. Ecker, J. Gasser, H. Leutwyler, A. Pich, E. de Rafael, Phys. Lett. B223 (1989) 425. 

K. Kampf, J. Novotny, J. Trnka, Eur. Phys. J. C50 (2007) 385-403. [ hep-ph/0608051] . 

K. Kampf, J. Novotny, J. Trnka, Phys. Rev. D81 (2010) 116004. [arXiv:0912.5289l 
[hep-ph]]. 

E. Witten, Nucl. Phys. B 223 (1983) 422. 

J.A. Schouten, Proc. Kon. Ned. Akad. v. Wet. 41 (1938) 709-716 



A. Nyffeler, Phys. Rev. D79 (2009) 073012. |arXiv:0901. 11721 [hep-ph]]; A. Nyffeler, PoS 
CD09 (2009) 080. [ arXiv:0912.1441l [hep-ph]]. 

G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22 (1980) 2157. 

S. J. Brodsky and G. P. Lepage, Phys. Rev. D 24 (1981) 1808. 

B. Aubert et al. [The BABAR Collaboration], Phys. Rev. D 80 (2009) 052002 
|arXiv:0905.4778l [hep-ex]]. 

S. S. Agaev, V. M. Braun, N. Offen, F. A. Porkert, Phys. Rev. D83 (2011) 054020. 
[arXiv:1012.467T1 [hep-ph]]. 

A. E. Dorokhov, [ arXiv:1003.4693l [hep-ph]]. 



K. Kampf and B. Moussallam, Eur. Phys. J. C 47 (2006) 723 |arXiv:hep-ph/0604125|. 



I. Larin et al. [PrimEx Collaboration], arXiv: 1009. 1681 [nucl-ex]. 

K. Abe et al. [ Belle Collaboration ], [ hep-ex/0610022] . 

K. Nakamura et al. [Particle Data Group], J. Phys. G 37 (2010) 075021. 

J. F. Donoghue, B. R. Holstein and Y. C. R. Lin, Phys. Rev. Lett. 55 (1985) 2766 
[Erratum-ibid. 61 (1988) 1527]. 

J. Bijnens, A. Bramon and F. Cornet, Phys. Rev. Lett. 61 (1988) 1453. 

J. Bijnens and K. Kampf, Nucl. Phys. Proc. Suppl. 207-208 (2010) 220 [arXiv:1009.5493l 
[hep-ph]]. 

K. Kampf, M. Knecht, J. Novotny and M. Zdrahal. larXiv:1103.098"2l [hep-ph]. 

F. Jegerlehner, A. Nyffeler, Phys. Rept. 477 (2009) 1-110. [arXiv:0902.3360l [hep-ph]]. 



41 



[52] J. Bijnens, E Pallante, J. Prades, Phys. Rev. Lett. 75 (1995) 1447-1450. 
[hep-ph/9505"25l]; J. B ijnens, E. Pallante, J. Prades, Nucl. Phys. B474 (1996) 379- 
420. |hep-ph/9511388| ; M. Knecht, A. Nyffeler, Phys. Rev. D65 (2002) 073034. 
[hep-ph/0111058] . 

A. Bertin et al. [ OBELIX Collaboration ], Phys. Lett. B414 (1997) 220-228. 

T. Gherghetta, J. I. Kapusta, T. M. Kelley, Phys. Rev. D79 (2009) 076003. 
|arXiv:0902.1998l [hep-ph]] . 

L. Cappiello, O. Cata, G. D'Ambrosio, |arXiv:1009.fT6T1 [hep-ph]]. 
S. Z. Jiang and Q. Wang, Phys. Rev. D 81 (2010) 094037 |arXiv:1001.03T5l [hep-ph]]. 
R. Unterdorfer and H. Pichl, Eur. Phys. J. C 55 (2008) 273 [arXiv:0801.2482 [hep-ph]]. 
O. Strandberg, [hep-ph/0302064| . 

A. A. Poblaguev et al, Phys. Rev. Lett. 89 (2002) 061803 [arXiv:hep-ex/0204006l . 

R. Kaiser, H. Leutwyler, Eur. Phys. J. C17 (2000) 623-649. [ hep-ph/0007l"0l] . 

G. W. Bennett et al. [ Muon G-2 Collaboration ], Phys. Rev. D73 (2006) 072003. 
[hep-ex/0602035] . 



[62] see http://www.g-2.bnl.gov/ 



42 



